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In this paper, we look at the following question. We consider cellular automata in the hyperbolic 
plane, see ||5] [16] [7] [9) and we consider the global function defined on all possible configurations. Is 
the injectivity of this function undecidable? The problem was answered positively in the case of the 
Euclidean plane by Jarkko Kari, in 1994, see Q. In the present paper, we show that the answer is 
also positive for the hyperbolic plane: the problem is undecidable. 

1 Introduction 

The global function of a cellular automaton A is defined in the set of all configurations. Note that when 
we implement an algorithm to solve a given problem, the initial configuration is usually finite. The study 
of the global function starts from another point of view. 

In the case of the Euclidean plane, the definition of the set of configurations is very easy: it is Q , 
where Q is the set of states of the automaton. 

In the hyperbolic plane, see J9l|7l, we have the following situation: we consider that the grid is the 
pentagrid or the ternary heptagrid, see (9). We fix a tile, which will be called the central cell and, around 
it, we dispatch a sectors, a G {5,7}: a = 5 in the case of the pentagrid, a = 7 in the case of the ternary 
heptagrid. We assume that the sectors and the central cell cover the plane and the sectors do not overlap, 
neither the central cell, nor other sectors: call them the basic sectors. Denote by J£« the set constituted 
by the central cell and a Fibonacci trees, see [9] [TH > each one spanning a basic sector. Then, a 
configuration of a cellular automaton A in the hyperbolic plane can be represented as an element of Q^ a , 
where Q is the set of states of A. If /a denotes the local transition function of A, its global transition 
function Ga is defined by: Ga{c){x) = f{c(x)), where c runs over Q^ a and x G & a . 

The injectivity problem for a cellular automaton consists in asking whether there is an algorithm 
which, applied to a description of /a would indicate whether Ga is injective or not. 

The goal of this paper is to prove the following result: 

Theorem 1 There is no algorithm to decide whether the global transition function of a cellular automa- 
ton on the ternary heptagrid is injective or not. 

It gives a negative answer to the previous question. Although this answer is not surprising, it is not 
at all obvious to prove it. As in the Euclidean case, the proof relies on the undecidability of the tiling 
problem. This is a well known result for the Euclidean plane, established by R. Berger in 1966, see HI, 
this proof being much simplified by R. Robinson in 1971, see iTTTl . The same result for the hyperbolic 
plane, this question was raised by R. Robinson in his 1971 paper was recently only established by the 
present author and J. Kari, in 2007, independently of each other, the proofs being completely different. 
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Figure 1 : The illustration of the new setting of the construction leading to the proof of the undecidability 
of the tiling problem in the hyperbolic plane. 

In this paper, we give a short account of the proof of the undecidability of the tiling problem, follow- 
ing the proof I established in lfT0l[T3l l8ll6l. in a new setting I indicated in lfT2l[TTi . This will be Section 2 
of this paper. Section 3 will introduce the mauve triangles, a new object on which the path, defined in 
Section 4, can be constructed, yielding the proof given in Section 5. We shall shortly discuss the situation 
opened by this result in our conclusion. 

We shall not repeat the construction of the interwoven triangles on which the construction of the 
mauve triangles relies. In section 2, we remind the basic properties of the mauve triangles and we 
introduce new ones. In section 3, we more carefully describe the construction of the path based on the 
mauve triangles. In section 4, we show how to derive the proof of theorem Q] 

2 The undecidability of the tiling problem 

In this section, we just introduce the new setting and then, sketchily remember the guidelines of the proof 
which is explained in |[T0l[T3l[T2l . 

The new setting also takes place on the ternary heptagrid, the tiling {7,3} of the hyperbolic plane: it 

In 

is based on the replication of a regular heptagon with — as the interior angle by reflection in its sides 
and then, recursively, by reflection of the images in their sides. 

We introduce four colours for the tiles: green, yellow, blue and orange which will be denoted by G, 
Y, B and O respectively. We attach the following substitution rules, see ITT21 ITTTl . as illustrated by the 
left-hand side picture of Figure [T] 

G -> YBG, Y -> YBG, B^BO,0^ YBO. 

As noted in ifTTTl . this allows to define two types of Fibonacci trees, see ||5] for the definition. If 
we define B to be a black node and G, Y and O to be white nodes, we obtain a central Fibonacci tree, 
see HJ[9l. If we define G and Y to be black and B with O to be white, then we have standard Fibonacci 
trees. We define seeds as G-nodes on an even isocline whose fatheris a 7-node and whose grandfather is 
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Figure 2: Construction of the interwoven triangles in the Euclidean plane: look at the green signal. 

a G-node. In Figure \T\ the tree rooted at a seed, in red in the figure, implements the representation of the 
trees of the mantilla, see |[T0l[T3l[T2l . We use this latter interpretation which allows us to construct the 
isoclines which are the horizontals of our constructions, both for the undecidability of the tiling problem 
and for the proof of TheoremQ] This is illustrated by the right-hand side picture of FigureQ] The convex 
arcs are attached to the black nodes, G and Y, and the concave ones are attached to the white nodes, here 
B and O. 

In this new setting, the isoclines are periodically numbered from to 3. We conclude this section by 
reminding the principle of the interwoven triangles, which is illustrated by Figure |2] The generation 
is in light blue in the figure. The even generations are in dark blue. Triangles, thick strokes, alternate 
with phantoms, thin strokes, for each generation. Triangles and only them generate the vertices and the 
bases of the triangles and phantoms of the next generation: light and dark blue triangles generate the red 
generations. The red triangles generate the dark blue generations. The reader is referred to iflOl [T3l [H HI 
for more details on this construction and its implementations in the hyperbolic plane. 

3 The mauve triangles 

The mauve triangles are constructed upon the interwoven triangles. 

The construction of the interwoven triangles needs a lot of signals, which entails a huge number of 
tiles, still around 9,000 of them in this new context, not taking into account the specific tiles devoted to 
the simulation of a Turing machine. The construction of this paper requires much more tiles, but we shall 
not try to count them. 

The mauve triangles are constructed from the red triangles of the interwoven ones. Any mauve 
triangle is triggered by a red triangle and conversely. The vertex of a mauve triangle T is that of a red 
triangle R. Its legs follow those of R. After the corner of R, they go on, on the same ray, until they meet 
the basis of the red phantoms which are generated by the basis of R. At this meeting, the legs meet the 
basis of the mauve triangle which coincide with the basis of the just mentioned red phantoms. In (H, we 
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Figure 3: An illustration of the mauve triangles. Inside the big triangle, from left to right, the 0-, 1-, 2- 
and 3-triangles of the previous generation. 

thoroughly describe that this construction can be forced by a finitely generated tiling. We refer the reader 
to these papers. 

3.1 Intersection properties of the mauve triangles 

As it stems from a red triangle, we say that a mauve triangle of the generation n is constructed on a red 
triangle of the generation 2n+l. Later, it will be useful to recognize the mauve triangles of generation 0: 
they have the colour mauve-0 and will be called mauve-0 triangles. 

From the doubling of the height with respect to the red triangles, the mauve triangles lose the nice 
property that the red triangles are either embedded or disjoint. This is no more the case for the mauve 
triangles. However, the overlappings and intersections of mauve triangles can precisely be described. 

The figure indicates the first points, FP, the mid-points, MP, and the low points, LP, of a mauve 
triangle T of generation n+l. On Figure |3l we can see that the 3-triangles cut the basis of T at their LP's 
and that the LP's of the 2-triangles are aligned with those of T. The isoclines which pass through the 
LP's, MP's and LP's of T cut smaller triangles at their LP's and for triangles of generation n— 1 if any, 
they are cut by these lines if and only if they are 2-triangles. The 3-triangles of the generation n which 
contains the vertex of T , if any, is called the hat of T . Note that inside T , its 0-triangles are all hatted. 

We can see that the intersection between mauve triangles are those which are indicated by Figure [3] 
and only those. A mauve triangle T intersects 3-triangles of the previous generation: always one of 
them at its basis, possibly another one near its vertex, at its legs. Also, T may be cut by the basis of a 
mauve triangle of a bigger generation, always at its LP's. We refer the reader to |[T4l for a proof of these 
properties. 

3.2 New notions: the /3-clines and the j8-, y-points 

From the fact that the 3-triangles 7}'s of the previous generation which are inside a mauve triangle T cut 
T along its basis and at their LP's, we can see that we can repeat this observation with the 3-triangles of 
the Tj's until we reach mauve-0 triangles. We call /3-cline the isocline of the mauve-0 triangles reached 
by this process and we say that the j3-cline is attached to T. Note that the same isocline is the /3-cline 
of all the 3-triangles which occur in the just considered construction. By definition, we call j3-points 
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of T the intersection of the legs of T with the /3-cline of the 2-triangles of the previous generation which 
it contains. Similarly, we call y-points of T the intersection of the legs of T with the /3-cline of their 
hat, if any. Note that if the hat of T does not exist, there are infinitely copies of T within the same 
set of isoclines which cross T which are hatted, so that for these mauve triangles the /3-cline exist and 
the considered isocline, the same for all these copies of T also intersect T, which allows to define the 
/-points of T. 

In fact, the y-points of T can be constructed in another way, thanks to the following remark. Assume 
that the hat H of T exists and assume that T belongs to the generation n+l. Consider a O-triangle Z of 
the generation n which is inside T. It is plain that the hat of Z exists if n > 0, as it is inside T and that it 
is a 3-triangle of the generation n— 1 which cuts the basis of H. Accordingly, H and the hat of Z have the 
same /3-cline. Consequently, to find the y-point of T, it is enough to take the intersection of the legs of T 
with the /3-cline of the hat of Z. Note that if Z is a mauve-0 triangle, the considered /3-cline is the isocline 
which passes through the vertex of Z. Consequently, this defines the y-points of T, independently from 
the existence or not of the hat of T. 

Note that the /3-, and y-points are not defined for mauve-0 triangles. 

All the points which we defined in Subsection 3.1 and in this one can be constructed by a small set 
of signals in terms of tiles. 

The constructions relies on the fact that the isoclines joining the FP's and the MP's as well as the 
basis of T have the following property, assuming that T belongs to the generation n+l: if the isocline 
meets a triangle of the generation i with i < n, it is a 2-triangle. This allows to locate the 0-, 1- and 
3-triangles of the generation n whose vertex is inside T. From this detection, we also can detect the 
2-triangles of the generation n which are inside T. 

Now, we have to indicate how to construct the LP's of T as well as its /3-cline. 

The construction of the LP's of T start from the MP's which send two signals: one along the leg of T 
and another along the isocline joining the MP's, in the direction of the other MP. When the horizontal 
signal meets the vertex of a red phantom, it goes down along it, until it reaches the mid-point M of the 
leg of the phantom. There, it goes back to the leg of T by following the isocline of M, circumventing all 
triangles which are met from their right-, left-hand side leg, the signal emanating from the left-, right- 
hand side leg of T. In this way, the first leg met by the signal whose laterality is that of the leg from 
which it came is the expected leg and it also meets the second signal sent by the MP. 

The construction of the /3-cline of T follows from this construction. A signal leaves each LP of T. 
They go to the opposite side. They follow the leg of T, then its basis and then, each time they meet a leg 
of a 3-triangle at its LP, they go down along this leg and, finding the corner, they go along the isocline 
of the basis of the corner, in the direction of the other corner. This signal is permitted only to go down 
or to follow an isocline. Consequently, the two considered signals will meet on the isocline which is the 
/3-cline of T, which provides the construction. 

The detection of a /3-cline can now be used to construct the /3- and y-points, combining the detection 
of a /3-cline with that of a 0-, or a 2-triangle. The pictures of Figure 0] illustrates the constructions of the 
/3- and the y-points of T. 

In |[T4l . we completely prove the following statement: 

Lemma 1 The mauve triangles together with the determination of their LP 's and mid-points can be 
constructed from a finite set of prototiles. 

which can be established using the above ideas. 
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Figure 4: Illustration of the construction of the jS- and y-points of a mauve triangle of a positive gen- 
eration. On the left-hand side: construction of the j3 -point. On the right-hand side, construction of the 
y-point. 

The y-point allows us to define a last couple of points on the legs of a mauve triangles, the HP's, i.e. 
high points. 

The HP's are defined with respect to the y-points as follows. If the /3-cline which crosses the y-point 
is a j3-cline of type 2, then the HP is the y-point. If this is not the case, then the HP of T is defined by 
the first intersection with a basis of a mauve triangle with the legs of T, when starting from the FP's and 
going up towards the vertex. If no basis is met before reaching the vertex, the HP of T is the intersection 
of the legs with the isocline which is just below that of the vertex of T. This is the case for the mauve-0 
triangles. 

4 A half-plane filling path 

Now, we turn to the construction of the path. 

As we shall see, the path has no cycle and most often, it consists of a single component which 
is plane-filling path. In one exceptional situation, the path breaks down into infinitely many infinite 
components. There is no cycle and each component fills up at least a half -plane. 

The points which we defined in the previous section can be seen as signboards which are placed on 
the path in order to indicate it in which regions it passes and which direction it has to follow. 

Our first task consists in describing these regions. In a second step, we shall see how the path fills up 
the regions, on the basis of an inductive construction. 

We define two types of basic regions. The first one is the set of tiles defined by a mauve triangle: 
its borders and its inside. Remember that the basis of a mauve triangle contains more than the majority 
of tiles resulting from the just given definition. It is considered as a basic region as once the path enters 
a mauve triangle T, it fills up T almost completely before leaving T. In fact, there is a restriction and 
the path fills a bigger area. Remember that the basis of T is crossed by the legs of a lot of triangles, 
the 3-triangles of the previous generation whose vertices are inside T and a lot of 2-triangles of smaller 
generations. The path enters each such triangle M through one of its LP's and it fills it up completely, 
leaving it through the HP of M which lies on the other leg of M. After a while, as will be seen, the path 
goes back to the basis of T and will meet the next small triangle which crosses this basis. Of course, 
inside a small triangle filled up by the path, the process is recursively repeated: if the basis of M is 
crossed by legs of smaller triangles, the path fills them up before going on on the basis of M. 
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Figure 5: A schematic representation of the path: 

On the left-hand side, inside a mauve-0 triangle. On the right-hand side, in between two consecutive 
mauve-0 triangles within the same latitude when the vertices belong to the same basis. 



This process allows us to define the latitude of a mauve triangle T. For a mauve-0 triangle, the lower 
border is the isocline of the basis of T and it belongs to the latitude. The upper border is the isocline 
of the vertex of T and it belongs to another latitude. For higher generations, the border of the latitude 
belonging to the latitude is defined by the isocline of the basis of T with this exception that each time 
this isocline is crossed by legs of smaller triangles, the border goes down along the legs and then follows 
the basis of the smaller triangles, repeating this way recursively, until a mauve-0 triangle is reached. 

In Figure|5] the picture (a) illustrates the path inside a mauve-0 triangle and the picture (b) illustrates 
the path in between two consecutive mauve-0 triangles inside the same latitude. Now, for this situation 
in between two consecutive mauve triangles of the same latitude, the path fills up the space between the 
upper and the lower borders. The situation is in general different from the picture (b). Figure [6]describes 
the general situation for consecutive mauve triangles of generation 1 standing within the same latitude. 
There, we can see the various configurations which may occur in between consecutive mauve-0 triangles 
of the same latitude. 

In Figure |6]we make use of the schematic representation used the pictures (c) and (d) of Figure[5]to 
represent the pictures (a) and (b). This is repeated in Figure [7J which represents the path inside a mauve 
triangle of generation 1 . 

Other details attached to the path are dealt with in lfl4l and we refer the reader to this paper. 

Now, we notice that Figures [6] and [7J can also illustrate the passage from the generation n to the 
generation n+l. We draw the attention of the reader on the following point. Due to the position of the 
entry and the exit points for the path in a triangle, the path fills up the bottom part of the triangle which 
is below the isocline joining its LP's in an almost cycle and then it crosses the other three strips roughly 
corresponding to 2-, 1- and O-triangles of the previous generation, in changing its direction each time it 
enters a new strip. We have a similar feature in between two consecutive triangles of the same latitude. 
This dictates a difference of behaviour when going to a strip to another. Now, this is regulated by the 
jS-clines of type 2: they delimit the bottom part of inside a triangle and they delimit the upper one in 
between two consecutive triangles of the same latitude. As by definition, the intersection of a /3-cline of 
type 2 with the legs of a triangle define the HP of this triangle, this allows to give a passage to the path 
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Figure 6: The path in between two triangles of generation 1. 




Figure 7: The path inside a triangle of generation 1. 

On the left-hand side: a 0- or a 1-triangle. On the right-hand side, a 3-triangle or a 2-triangle which is 
cut by a mauve triangle of a bigger generation. 
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for going back to the origin of the almost cycle. The intersection of a /3-cline which is not of type 2 with 
the legs of a triangle prevent the path to enter the triangle or to exit from it. 

With these indications, we conclude that the construction is correct and that the path possesses the 
following property: 

Lemma 2 The path contains no cycle. 

Lemma 3 For any tile X, the path on one side of z fills up infinitely many mauve triangles of increasing 
sizes. 

From these lemmas, as established in lfl4l . we get: 

Corollary 1 If there are only finite basic regions, the path goes through any tile of the plane. 

Now, there is an exceptional situation in which we have infinitely many infinite triangles. In fact, as 
the path fills up all triangles, the infinite ones included, it splits into infinitely many components which 
still have no cycle and which fill up the region contained in an infinite triangle and the one in between 
this infinite triangle and the next one in the same infinite latitude. Now, for these components, Lemmas [2] 
and [3] are still true, see lPT4l . 

5 Proof of the main theorem 

We can now prove Theorem Q] 

The proof follows the argument of Q, with an adaptation to the case of infinite triangles. 

We define a direction for the path constructed in Section 4. To this aim, we introduce three hues in 
the colour used for the signal of the path. One colour calls the next one and the last one calls the first 
one. The periodic repetition of this pattern together with the order of the colours define the direction. 
This allows to define the successor of a tile on the path, which we formalize by a function 8 from 2Z to 
the tiling such that 8(n+l) is the successor of 8{n) on the path. 

Consider M a deterministic Turing machine with a single head and a single bi-infinite tape which is 
assumed to be initially empty. From lITOl IT3~ll . we can define a finite set of tiles Tm such that Tm tiles the 
hyperbolic plane if and only if M does not halt. When Tm tiles the hyperbolic plane, it defines infinitely 
many triangles of infinitely many sizes and in each triangle, the simulation of M is resumed from its 
beginning, running until the basis of the triangle is reached. An automaton Am is attached to M and its 
states are defined by D x {0, 1} x Tm, where D is the set of tiles which defines the tiling which we have 
constructed in Sections 3 and 4. The 0, 1 -component of a state is called its bit. We can still tile the 
plane as the tiles of Tm are ternary heptagons but the abutting conditions may be not observed: if it is 
observed with all the neighbours of the cell x, the corresponding configuration is said to be correct at x, 
otherwise it is said incorrect. When the considered configuration is correct at every tile for D or at every 
tile for Tm, it is called a realization of the corresponding tiling. 

As in O, the transition function does not change neither the D- nor the T^-component of the state of 
a cell x: it only changes its bit. As in |3J, we define Am(c(x)) = c(x) if the configuration in D or in T is 
incorrect at the considered tile. If both are correct, we define Am{c{x)) = xor (c(x),c (8 (x))). It is plain 
that if M does not halt, Tm tiles the hyperbolic plane and there is a configuration of D and one of Tm 
which are realizations of the respective tilings. Then, the transition function computes the xor of the bit 
of a cell and its successor on the path. Hence, defining all cells with and then all cells with 1 define two 
configurations which Am transform to the same image: the configuration where all cells have the bit 0. 
Accordingly, Am is not injective. 
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Conversely, if A M is not injective, we have two different configurations cq and c\ for which the image 
is the same. Hence, there is a cell x at which the configurations differ. Hence, the xor was applied, which 
means that D and T are both correct at this cell in these configurations and it is not difficult to see that the 
value for each configuration at the successor of x on the path must also be different. And so, following 
the path in one direction, we have a correct tiling for both D and Tm- Now, from Lemma |3l as the path 
fills up infinitely many triangles of increasing sizes, this means that the tiling realized for Tm is correct in 
these triangles. But, in each triangle, the computation of M is simulated from its beginning. Accordingly, 
the Turing machine M never halts. And so, we proved that Am is not injective if and only if M does not 
halt. Accordingly, the injectivity of Am is undecidable. | 

6 Conclusion 

In the Euclidean plane, the surjectivity of the global function of cellular automata is tightly connected 
with the injectivity of the function. This is not the case in the hyperbolic plane, see |[T5l . so that the 
question of the surjectivity is open in this case. Note that the present construction can be generalized to 
any grid {p,3} of the hyperbolic plane, with p > 7. 
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